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ABSTRACT 

We present a comparison of histograms of CO (2-1) line centroid velocity increments 
in the p Ophiuchi molecular cloud with those computed for spectra synthesized from 
a three-dimensional, compressible, but non-starforming and non-gravitating hydrody- 
namic simulation. Histograms of centroid velocity increments in the p Ophiuchi cloud 
show clearly non-Gaussian wings, similar to those found in histograms of velocity in- 
crements and derivatives in experimental studies of laboratory and atmospheric flows, 
as well as numerical simulations of turbulence. The magnitude of these wings increases 
monotonically with decreasing separation, down to the angular resolution of the data. 
This behavior is consistent with that found in the phase of the simulation which has 
most of the properties of incompressible turbulence. The time evolution of the mag- 
nitude of the non-Gaussian wings in the histograms of centroid velocity increments in 
the simulation is consistent with the evolution of the vorticity in the flow. However, we 
cannot exclude the possibility that the wings are associated with the shock interaction 
regions. Moreover, in an active starforming region like the p Ophiuchi cloud, the effects 
of shocks may be more important than in the simulation. However, being able to identify 
shock interaction regions in the interstellar medium is also important, since numerical 
simulations show that vorticity is generated in shock interactions. 

Subject headings: hydrodynamics - ISM: clouds - ISM: kinematics and dynamics - 
turbulence 
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1. Introduction 

Early spectroscopic observations of interstellar lines 
of Hi, OH, and CO have revealed that observed line 
widths (or velocity dispersions) in interstellar clouds 
are larger them thermal line widths expected for these 
low-temperature regions (see e.g. Myers 1997 and 
references therein). These large line widths are in- 
dicative of supersonic motions of the gas, however, 
the exact nature of these motions is still a subject 
of controversy. Proposed explanations generally in- 
volve turbulent motions of the gas, be it hydrody- 
namic or magneto-hydrodynamic (e.g. Scalo 1987; 
Falgarone 1997). 

The idea that the gas motions in interstellar clouds 
are of turbulent nature was first proposed by Zuck- 
erman & Evans (1974). As discussed by Falgar- 
one (1997), gas motions in interstellar clouds exhibit 
many features of incompressible turbulence. One of 
them is the size — line-width relation first discussed by 
Larson (1981). This power-law scaling with an expo- 
nent of ~0.3-0.5j found in the velocity field of molec- 
ular clouds is reminiscent of the Kolmogorov cascade 
in incompressible turbulence. Another possible sign- 
post of turbulence is the fractal structure of the in- 
tensity contours in molecular clouds (e.g. Scalo 1990; 
Falgarone, Phillips, &: Walker 1991). In both cases 
alternative explanations have been proposed. An ex- 
planation for the size — line-width relation observed 
from cloud to cloud in a given molecular line has been 
proposed as based on virial equilibrium and the rela- 
tively small variation of the observed column density 
between clouds of different size (Myers 1997). Expla- 
nations for the observed scaling laws have also been 
given in terms of the framework of magnetically sup- 
ported cloud models (Scalo 1987). Alternative ex- 
planations for the fractal cloud geometry based on 
collisional fragmentation or gravitational clustering 
and virial balance have also been proposed (Nozakura 
1990; Pfenniger & Combes 1994). 

One of the tools employed in studies of gas motions 
in the interstellar medium is the analysis of shapes of 
molecular line profiles. Falgarone & Phillips (1990) 
argued that the non-Gaussian CO line wings in inac- 
tive regions without associated star formation activity 
represent a direct observational signature of the tur- 
bulent nature of the gas flow within molecular clouds 
and of the existence of regions of intermittent tur- 
bulent activity. In a subsequent study Falgarone et 
al. (1994) showed that line profiles synthesized from a 


three-dimensional turbulent, compressible, but non- 
starforming and non- gravitating simulation (Porter, 
Pouquet & Woodward 1994) are in fact statistically 
similar to the CO line profiles observed in quiescent 
molecular clouds. However, Dubinski, Narayan, &: 
Phillips (1994) showed that non-Gaussian line pro- 
files can also be produced from a random velocity 
field with a Kolmogorov power spectrum. 

It is known from numerical simulations and atmo- 
spheric and laboratory measurements (e.g. Anselmet 
et al. 1984; Gagne 1987; Vincent & Meneguzzi 1991) 
that intermittency of turbulence manifests itself through 
non-Gaussian wings in probability density functions 
(PDFs or histograms) of velocity increments and deriva- 
tives . These kinds of measurements are not directly 
possible in the case of the interstellar medium where 
one has information integrated over a line-of-sight col- 
umn defined by the passage of the telescope beam 
through the medium. However, in a recent study (Lis 
et al. 1996), motivated by observational work of Mi- 
esch & Scalo (1995), we showed that histograms of 
the centroid velocity increments (or differences be- 
tween line centroid velocities at positions separated 
by a given distance), for sections of the simulation 
data cube corresponding to the ISM columns, also 
show non-Gaussian wings, similar to those found in 
experiments and numerical simulations of incompress- 
ible turbulence. Due to line-of-sight averaging, the 
wings seen in PDFs of centroid velocity increments 
are not as pronounced as those in PDFs of veloc- 
ity increments calculated over the w r hole data cube. 
Nevertheless, the effect is clearly present. We also 
demonstrated that the lines of sight contributing to 
the non-Gaussian wings of the PDFs of centroid ve- 
locity increments trace a filamentary structure, which 
follows the distribution of the two vorticity compo- 
nents involving cross-derivatives of the line-of-sight 
component of the velocity field. This suggests that 
the wings are a manifestation of the turbulent nature 
of the flow. 

The method presented by Lis et al. (1996) provides 
a useful new tool for studying the effects of intermit- 
tency of turbulence on physics and chemistry of the 
interstellar medium. In the present paper wo apply 
this method to a large-scale map of the CO (2-1) 
emission from the p Ophiuchi molecular cloud. It is 
true that this region with active low-mass star forma- 
tion may not be ideal for a study of turbulence, since 
the shocks associated with embedded young stellar 
objects can have their signatures in the velocity field. 
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On the other hand, the energy injected into the cloud 
by the embedded sources may also play an important 
role in generating the turbulent cascade. A statisti- 
cal comparison of the velocity field in both active and 
quiescent regions is thus required. (Observations of a 
truly quiescent region are in progress.) In addition, 
we extend our analysis to two additional epochs of the 
simulation of Porter et al. (1994) to further investi- 
gate the origin of the non-Gaussian wings in the PDFs 
of centroid velocity increments in the simulation. 

2. Observations 

Observations of the CO (2—1) emission from the 
p Oph molecular cloud presented here were carried 
out between 1994 July and 1995 March using the 
10.4-meter Leighton telescope of the Caltech Submil- 
limeter Observatory (CSO) on Mauna Kea, Hawaii. 
The data set consists of 18 partially overlapping on- 
the-fly (OTF) maps. The relative calibration between 
the various maps was determined from the overlap re- 
gions and the maps were scaled by the corresponding 
intensity ratios. The final map thus has a uniform 
calibration. In the OTF technique not all the spec- 
tra are taken on the same regular-grid in the equato- 
rial coordinate system, a final regular grid map with 
15" spacing was produced by averaging spectra falling 
with 8" radius from a given grid point [the CSO beam 
at the frequency of the CO (2-1) line is ~31"]. This 
final resampled map consists of 8970 spectra taken 
over an area of ~30 x 27'. A gray-scale image of the 
integrated line intensity between —4 and 12 kms 
is shown in Figure 1 and velocity channel maps are 
shown in Figure 2. In spite of the fact that the 
data were generally taken under less than optimum 
weather conditions, the resulting signal-to-noise ratio 
for the integrated intensity is quite high (12-280). In 
the subsequent analysis we use a 15" pixel as a length 
unit for the x and y axes. The size of our map is 122 
by 109 pixels. 

3. Discussion 

The molecular cloud near p Oph is one of the 
closest clouds, ranging in distance from about 80 to 
200 pc. This thickness is thought to be real (de Geus, 
de Zeeuw & Lub 1989; de Geus k Burton 1991). The 
cloud is brightly illuminated by the Sco OB 2 asso- 
ciation, the closest such association to the Sun. The 
center of the Upper-Scorpius subgroup lies to the west 
of the p Oph cloud and its stars lie at the same dis- 


tance or beyond the cloud. Some of the stars appear 
to lie within the cloud. 

In addition, the p Oph cloud is a site of ongoing 
low-mass star formation. This cloud has been the 
subject, over the years, of many near-infrared surveys 
(e.g., Grasdalen, Strom k Strom 1973; Elias 1978; 
Wilking & Lada 1983; Greene k Young 1992; Bar- 
sony et al. 1997). Most of these surveys have been 
of limited sensitivity or of small areas. In order to 
decide on the membership of a source in the embed- 
ded infrared stellar cluster, it is good to have data 
in a variety of wavelength bands, optical, near-, mid- 
, far-infrared, radio, and X-ray. Only of order 100 
sources have been identified so far as being probable 
p Oph cloud members (Wilking et al. 1983; Strom, 
Kepner k Strom 1995; Casanova et al. 1995). The 
rest simply do not have enough information available 
for them. The most recent near-infrared survey, Bar- 
sony et al. (1997), is a high sensitivity survey of a 
large area (~ 1° square), which catalogs 4495 sources. 
Analysis of this data set in terms of cloud member- 
ship has not been completed. Probable p Oph cloud 
members are identified in Figure 3. The symbols in- 
dicate near-infrared, X-ray and radio sources and op- 
tical emission line stars. Many sources have been de- 
tected in more than one band. The sources are seen 
to be concentrated in two regions - near the dense 
molecular cores, called p Oph A and C. This is partly 
due to the fact that these areas are among the most 
thoroughly surveyed, but nevertheless appears to be 
a real effect. 

The most luminous of the stars in the region 
that we have mapped are HD 147889 (B2 V; 5500 
L 0 ; 16 h 22 m 22?8, -24°21'07"), SI (B3 V; 1500 L©; 
16 h 23 m 32?8, -24°16'44") and SR 3 (B9 V; 125 L©; 
16 h 23 m 07!7, -24°27'26"). There are several sources 
with luminosities in the 10-50 L© range but most of 
them have luminosities of order 1 L© or less (Wilking, 
Lada k Young 1983). 

Stars can input energy into molecular clouds in 
ways other than through optical radiation. Many of 
the embedded stars are X-ray sources (Montmerle et 
al. 1983; Casanova et al. 1995). Also, there are sev- 
eral outflow sources known. Among these are the 
source of the extremely long outflow seen in Fig- 
ure 1 - VLA 1623 (16 h 23 m 25?0, -24°17'47"; Andre, 
Ward-Thompson, k Barsony 1990). Also, SR 4 
(16 h 22 m 54?8, -24°14'01") has been found to be the 
source of a string of HH objects (Wilking et al. 1997). 
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3.1. Centroid Velocity Increments in the p Ophi- 
uchi Cloud 

A gray-scale map of the CO (2-1) centroid veloc- 
ity in the p Oph cloud is shown in the top panel of 
Figure 4. In the bottom panels we show PDFs of 
the centroid velocity increments computed for sepa- 
rations A = 2, 5, 10, and 20 pixels (1 pixel = 15" = 

0.5 x FWHM), As in Lis et al. (1996), all histograms 
presented here are normalized to a zero mean and a 
standard deviation of unity. A normalized Gaussian 
distribution is also shown as a dotted line in each 
panel. The magnitude of the non-Gaussian wings in 
the histograms of the centroid velocity increments de- 
creases monotonically with the increasing separation 
A. This evolution is similar to that in the simulation 
(see Fig. 6), although for the p Oph histograms, the 
wings are somewhat more pronounced at small sepa- 
rations and the core of the histogram is more peaked. 

Several morphological features can be identified in 
the centroid velocity map. The most prominent is 
a narrow low-velocity (blue-shifted) filament running 
at ~45° position angle in the north-east quarter of 
the map. Emission associated with this filament can 
be seen at velocities between -1 and 1 Ions'" 1 in the 
channel maps (Fig. 2). This feature is the blue lobe 
of the molecular outflow originating near VLA 1623 
(Andre et al. 1990). A somewhat less prominent 
high-velocity (red-shifted) filament can be seen to the 
south-west of the low- velocity filament in the centroid 
velocity map (Fig. 4). No feature corresponding to 
the high-velocity filament can be distinguished easily 
in the integrated intensity map. In addition, two com- 
pact regions with low and high velocities are seen on 
the opposite sides of the two narrow filaments. There 
is also an east- west filament in the north-west quarter 
of the map. All these features can be identified in the 
channel maps (Fig. 2). 

3.2. Time Evolution of the Centroid Velocity 
Increments in the Simulation and the Ef- 
fects of Shocks 

As described in Falgarone et al. (1994), the time 
evolution of the simulation of Porter et al. (1994) is 
described as comprising three major phases. Their 
duration is conveniently expressed in terms of the 
acoustic time r ac = Lq/cs which is the sound cross- 
ing time of the energy containing scale, L 0 = /k 0 , 

where cs is the sound speed in the medium. The first 
short phase ends at t ~ 0.3r ac with the formation of 


first shocks. Shock interactions are the basic effects 
in the second supersonic phase where the vorticity 
generated by these interactions, as well as by vortex 
stretching, allows most of the kinetic energy to be 
transferred to solenoidal modes. The energy dissipa- 
tion rate is thus considerably reduced, which explains 
the long duration of this phase, up to t ~ 2.1r ac . In 
this phase and thereafter the flow has most of the 
properties of incompressible turbulence. The power 
spectrum for the solenoidal part of the flow (which 
contains more than 90% of the energy) is close to that 
predicted by Kolmogorov [£(£) ~ A;” 5 / 3 ]. The PDF 
of the vorticity in this phase shows a conspicuous non- 
Gaussian tail, illustrative of the intermittency of the 
flow. The last, post-supersonic, phase is that of the 
self-similar decay of the turbulent velocity spectrum. 

Lis et al. (1996) studied centroid velocity incre- 
ments in the second phase of the simulation, t = 
1.2r flC , when the solenoidal small-scale modes have 
already gained a large fraction of the energy and al- 
most all strong shocks have disappeared (Porter et 
al. 1994). To better understand the origin of the non- 
Gaussian wings in the histograms of centroid velocity 
increments, we have performed a similar analysis for 
the remaining epochs of the simulation. Results for 
the first phase, t = 0.5r oc , just after the shock forma- 
tion, are shown in Figure 5. At this early time, the 
vorticity is very low and small-scale disturbances have 
not yet developed. The velocity field is not yet char- 
acteristic of turbulence. The centroid velocity map 
for the second, supersonic, phase of the simulation 
(£ = 1.2r ac ) is shown in Figure 6 (top panel) along 
with the histograms of the centroid velocity incre- 
ments (bottom panels). Finally, results for the third 
phase of self-similar decay of turbulence (t = 2.4r ac ) 
are shown in Figure 7. 

In the first phase of the simulation (Fig. 5) non- 
Gaussian wings in the histograms of centroid veloc- 
ity increments are seen only for small separations, A. 
For A £ 15, the histograms have essentially Gaus- 
sian shapes. As time progresses, the magnitude of 
the wings increases. The effect is most pronounced 
at larger separations (A £ 7); non-Gaussian wings 
are clearly presents even for A = 30 in the second 
phase of the simulation (Fig. 6). The magnitude of 
the wings increases monotonically with decreasing A 
down to the numerical resolution of the simulation (~ 
several pixels). At even later times (Fig. 7), the mag- 
nitude of the wings increases further for large separa- 
tions (A ;> 15). However, the wings in the histograms 
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for small separations are not very prominent. As dis- 
cussed by Falgarone et al. (1994), in this late phase of 
the simulation most of the structure is at small scales, 
while large scale structures have already disappeared. 
The line-of-sight averaging may thus be more effective 
in smearing out the most extreme events for small sep- 
arations than at the earlier epochs of the simulation. 

The time evolution of the non-Gaussian wings in 
the histograms of centroid velocity increments in the 
simulation (in particular the increase in the magni- 
tude of the wings between phase one and two) is con- 
sistent with the evolution of vorticity in the flow, as 
described above. However, shocks may also have their 
signatures in the velocity field in the supersonic phase 
of the simulation. We can show based on a simple ge- 
ometrical argument that planar shocks will have lim- 
ited influence on the shapes of the histograms of cen- 
troid velocity increments. Lines of sight contributing 
to the non-Gaussian wings are associated with large 
centroid velocity shifts on small angular scales. In 
the case of a planar shock with the velocity in the 
plane of the sky (perpendicular to the line of sight) 
the shock front is well defined in the observations. 
However, there is no line-of-sight velocity jump, since 
the velocity of the shock is in the plane of the sky. Al- 
ternatively, for a shock with the velocity along the line 
of sight, there is a large line-of-sight velocity jump be- 
tween the pre-shock and post-shock gas, but we do not 
observe a well defined shock front (the pre-shock and 
post-shock gas overlap spatially on the sky). There- 
fore, in neither case will there be a strong contribu- 
tion to the centroid velocity increments. However, 
the regions of interactions between shocks are likely 
to exhibit large velocity gradients over small angu- 
lar scales. Therefore, we cannot exclude the possi- 
bility that shock interactions contribute to the non- 
Gaussian wings in the histograms of centroid velocity 
increments. However, being able to identify shock in- 
teraction regions is also interesting, since as discussed 
by Porter et al. (1994), the vorticity is generated in 
shock interactions. 

The increment method will also select boundaries 
of cloud fragments spatially superposed on the sky, 
which are of no interest to us and can only be identi- 
fied by inspecting individual spectra. 

3.3. Distribution of Points Selected by Incre- 
ments Method 

Lis et al. (1996) showed that the lines of sight con- 
tributing to the non-Gaussian wings of the PDFs of 


centroid velocity increments in the simulation trace 
a filamentary structure, which follows the distribu- 
tion of the two vorticity components involving cross- 
derivatives of the line-of-sight component of the ve- 
locity field. The distributions of the lines of sight 
in the p Oph map with centroid velocity increments 
greater than 3.5 for a separation A = 2 are shown 
in Figure 8. As in the simulation, the points are not 
randomly distributed, but instead trace a number of 
coherent structures described in § 3.1, namely the nar- 
row filaments, high- and low- velocity compact sources 
and the east-west filament. 

3.4. Optical Depth Effects 

The CO (2-1) lines in the p Ophiuchi cloud have 
high optical depths. In an optically thick case, the 
line emission does not trace uniformly the velocity 
field along the line of sight, but instead it is sensi- 
tive primarily to the gas near the surface of the cloud 
(note that the optical depth in the line wings may 
still be low). Unless the turbulent velocity field in 
the surface layer is statistically different from that 
inside the cloud, we do not expect significant differ- 
ences in the shape of histograms of centroid veloc- 
ity increments for optically thick and thin lines. Lis 
et al. (1996) considered centroid velocities calculated 
with density weighting (appropriate for optically thin 
lines), as well as uniform weighting (more appropri- 
ate for optically thick lines originating in a clumpy, 
supersonic medium) and showed that the resulting 
histograms are similar. 

To study the optical depth effects more quantita- 
tively, including spatial overlap of optically thick cells, 
we computed simulated spectra for the second epoch 
of the simulation for different optical depths through 
the cloud assuming LTE conditions. We assumed that 
the optical depth through each cell of the data cube is 
proportional to the gas density in the cell, and the ex- 
citation temperature is uniform over the whole cube. 
Figure 9 shows a comparison of histograms of cen- 
troid velocity increments computed for an increment 
A = 3 for the second phase of the simulation, in opti- 
cally thin and thick cases (average line center optical 
depths of 0.01 and 10, black and gray lines, respec- 
tively). As expected, the distributions corresponding 
to optically thin and thick cases are very similar, in- 
dicating that optical depth effects do not influence 
significantly the shapes of the histograms. From the 
point of view of statistical sampling of the molecu- 
lar cloud velocity field, using optically thick lines is 


5 



an advantage, since the emission is strong and it is 
relatively easy to collect a large enough data set for 
a meaningful statistical study. However, the distri- 
bution of points contributing to the wings of the his- 
tograms is different in the two cases. Since in the opti- 
cally thick case the emission does not trace uniformly 
the gas along the line of sight, the increment method 
selects preferentially the regions with increased vor- 
ticity located closest to the observer. High vorticity 
regions located along the line of sight deeper inside 
the cloud would largely remain undetected. There- 
fore, from the point of view of selecting lines of sight 
for comparative chemical and physical studies of tur- 
bulent regions with the bulk of the gas, the increment 
method is best applied to optically thin lines. 

4. Summary 

The time evolution of the non- Gaussian wings 
in the histograms of centroid velocity increments 
in a compressible, but non-star-forming and non- 
gravitating turbulent simulation is consistent with 
the evolution of the vorticity in the flow. In partic- 
ular, the magnitude of the wings increases between 
phase one and phase two of the simulation as shocks 
dissipate and the vorticity is generated. However, we 
cannot exclude the possibility that shock interactions 
contribute to the non-Gaussian wings present in the 
simulation. 

Histograms of CO (2-1) centroid velocity incre- 
ments in the p Ophiuchi cloud are similar to those 
found in the simulation. This is consistent with the 
proposition that ISM velocity structure is vorticity 
dominated like that of the turbulent simulation. How- 
ever, the region we have studied contains some ac- 
tive star formation, as indicated by the presence of 
infrared sources and molecular outflows. As a re- 
sult shocks may have more important effects on the 
velocity field structure and molecular line shap^ in 
this region compared to the simulation. In particu- 
lar, regions of shock interactions may contribute to 
the wings of the histograms of centroid velocity in- 
crements in the p Ophiuchi cloud. However, being 
able to identify shock interaction regions in the ISM 
is also important, since numerical simulations show 
that vorticity is generated in shock interactions. 

We have demonstrated that the increment method 
can be successfully applied as a tool for studying the 
velocity field in interstellar molecular clouds. The 
same analysis now has to be applied to a carefully se- 


lected quiescent region without signs of ongoing star 
formation. The presence of non-Gaussian wings in 
histograms of line centroid velocity increments in such 
a region would be a strong argument in favor of the 
turbulent nature of the velocity field in interstellar 
molecular clouds. Also, since the energy injected into 
the cloud by the embedded sources may play an im- 
portant role in generating the turbulent cascade, a 
comparison of the histograms of centroid velocity in- 
crements in active and quiescent regions may help in 
understanding similarities and differences in the sta- 
tistical properties of the velocity field between differ- 
ent regions in the ISM. 
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JPL Supercomputing Project which is sponsored by 
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Fig. 1. — Gray-scale map of the CO (2-1) integrated intensity between -4 and 12 kms -1 in the p Oph cloud. 
Contour levels are from 20 to 70 by 10 K kms -1 . The intensity has not been corrected for the main beam 
efficiency, which is ~75% at the frequency of the CO (2-1) line. 
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Fig. 2. — Channel maps of the CO (2-1) emission in the p Oph cloud. Center velocities are written in the upper-left 
corner of each panel. Contour levels are 5, 10, 15, 22, and 30 K. 




IR Sources 
X-ray Sources 
Radio Sources 
Hll Regions 
Emission Line Stars 



IR Sources 
X-ray Sources 
Radio Sources 
Hll Regions 
Emission Line Stars 


a 1t50 


Fig. 3. — Locations of near-infrared, X-ray, radio sources, and emission line stars associated with the p Oph cloud. 
Sources marked as Hll regions are bright compact radio sources from Brown & Zuckerman 1975. Remaining radio 
sources are mostly T-Tauri stars. Symbols are overlayed on top of images of the CO (2-1) integrated intensity in 
the velocity intervals 4-8 kins** 1 and 0-4 kms" 1 (upper and lower panels respectively). 


10 






-5 0 5 -5 0 5 

Normalized Centroid Velocity Increment 


Fig. 4. — (Top) Gray-scale map of the CO (2-1) line centroid velocity in the p Oph cloud. Spatial coordinates 
are in pixels (one pixel equals to 15", or one half of the FWHM beam-size). ( Bottom ) Normalized histograms of 
centroid velocity increments in the p Oph cloud for separations A = 2, 5, 10, and 20 pixels. 
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Fig. 5 — Map and histograms of the line centroid velocity increments for the first phase of the simulation (t — 
0.5r oc ). The non-Gaussian wings are less pronounced compared to phase two (Fig. 6) for any value of A. For 
A 15 histograms are close to a Gaussian. 
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Fig. 6. — (Top) Gray-scale map of the line centroid velocity for the second phase ( t = 1.2r ac ) of the simulation of 
Porter et al. (1994). Note that the axes are reversed compared to Fig. 6 of Lis et al. (1996). (Bottom) Normalized 
histograms of centroid velocity increments computed for separations A = 3, 7, 15, and 30 pixels. 
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Fig. 7. — Map and histograms of the line centroid velocity increments for the third phase of the simulation 
(£ = 2.4r oc ). The non-Gaussian wings are still present for large separations (A £ 15). For small separations 
histograms approach a Gaussian shape as small scales are dissipated in this phase of self-similar decay of turbulence. 
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pig g. — Distribution of the lines of sight in the Oph A centroid velocity map with normalized centroid velocity 
increments greater than 3.5 for a separation A = 2 pixels. 


5 


Centroid Velocity 



Normalized Centroid Velocity Increment 

Fig. 9. — Comparison of histograms of centroid velocity increments for optically thin and thick lines. The histograms 
were computed for a separation A = 3 for the second phase of the simulation. Black and gray lines correspond to 
a mean line center optical depths of 0.01 and 10, respectively. 
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